AD-A03J  609  COAST  6UAR0  RESEARCH  ANO  DEVELOPMENT  CENTER  6R0T0N  CONN  F/S  13/10 
HYDROSTATICS  ANO  EQUILIBRIUM  CHARACTERISTICS  COMPUTER  PROGRAM  F— ETC(U) 
APR  76  P SCCREST 

UNCLASSIFIED  C6R/DC-6/76  USC6-0-6S-76  NL 


• 

■W&r 

* 

|C. 

•- 

i,  - • '“*1 

mm 

rHRiini .{ 

I1  1 

_J 

i 

Ws::‘: 

' f'-: 

END 

DATE 

FILMED 

2-77 

5HDA033609 


••Mrt  No.  CG-D-65-76 


. 


HYDROSTATICS  AND  EQUILIBRIUM  CHARACTERISTICS 
COMPUTER  PROGRAM  FOR  RECREATIONAL  BOATS 


Peter  Secrest 


U.  S.  Coast  Guard  Research  and  Development  Center 
Avery  Point 

Groton,  Connecticut  06340 


Document  Is  available  to  the  public  through  the 
National  Technical  Information  Service, 
Springfield,  Virginia  22161 


Frapared  for 

EEMITNENT  IF  TRANSPORTATION 
•MTU  STATES  COAST  GOARS 

Off  I CO  of  Roaoorah  oad  Dovolopmont 
Woehlftfton,  D.C.  JOSOO 


OV 


NOTICE 

This  document  is  disseminated  under  the  sponsorship  of 
the  Department  of  Transportation  in  the  interest  of 
information  exchange.  The  United  States  Government 
assumes  no  liability  for  its  contents  or  use  thereof. 

The  contents  of  this  report  reflect  the  views  of 
the  U.S.  Coast  Guard  Research  and  Development  Center,  which 
is  responsible  for  the  facts  and  accuracy  of  data  presented. 
This  report  does  not  constitute  a standard,  specification 
or  regulation. 


i 


2 Government  Accession  No. 


Technical  Report  Documentation  Page 

3.  Recipient'*  Cotolog  No.  | 


HYDROSTATICS  AND  JRUILIBRIUM  CHARACTERISTICS 
COMPUTER  PROGRAM  FOR  RECREATIONAL  BOATS  m ^ 


Otmmg  Urgant  zd*iOn  Nome  and  Address 

U.  S.  Coast  Guard  Research  and  Development  Center 
Avery  Point 

Groton,  CT  06340  ^ 


12.  Sponsoring  Agency  Nome  ond  Address 

Department  of  Transportation 
U.  S.  Coast  Guard 

Office  of  Research  and  Development 
Washington,  DC  20590 


IS.  Supplementary  Notes 


CGR&DC  6/76 


16  Abstract 

is  report  is  intended  to  serve  as  complete  documentation  of  a versatile  computer 
program  which  calculates  hydrostatics  and  equilibrium  characteristics.  It  differs 
from  other  hydrostatics  programs  in  its  applicability  to  odd-shaped  boats  and  in 
its  capability  of  performing  equilibrium  calculations. 

The  integration  scheme  uses  trapezoidal  rule  throughout , allowing  random  spacing 
of  points  and  stations  defining  the  hull  form.  This  allows  almost  every  type  of 
hull  form  to  be  represented  and  computations  performed. 

The  equilibrium  characteristics  feature  allows  determination  of  transverse  and 
longitudinal  righting  moment  curves,  heave  restoring  forces,  trim  vs.  heel  angle 
coupling,  etc.  It  also  is  possible  to  determine  the  draft,  trim  and  heel  for 
any  loading  condition. 

Included  in  the  report  are  a program  explanation,  user  guide,  program  listings, 
program  flow  charts,  and  a sample  run. 
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1.0  PROGRAM  CHARACTERISTICS  OVERVIEW 


This  report  Is  Intended  to  serve  as  complete  documentation  of  a 
versatile  computer  program  which  calculates  hydrostatics  and  equilibrium 
characteristics.  Its  chief  differences  from  other  hydrostatics  programs 
are  applicability  to  odd-shaped  boats  and  capability  of  performing  equilibrium 
calculations. 

What  makes  this  program  applicable  to  odd-shaped  boats?  The  Integra- 
tion scheme  utilizes  trapezoidal  rule  throughout,  which  allows  the  points 
defining  the  hull  to  be  spaced  at  completely  random  intervals,  with  closer 
spacing  in  regions  of  greatest  curvature.  The  transverse  sections  through 
the  hull  can  also  be  spaced  at  random  intervals.  With  this  method  of 
hull  definition,  virtually  every  type  of  hull  form — Including  conventional 
ship  shape  hulls,  boats  with  hard  chines,  cathedral  hulls,  stepped  planing 
boats  and  catamarans — can  be  represented  and  computations  performed. 

The  equilibrium  characteristics  feature  allows  determination  of 
transverse  and  longitudinal  righting  moment  curves,  heave  restoring  forces, 
trim  vs.  heel  angle  coupling,  heave  vs.  heel  angle  coupling,  trim  vs. 
heave  coupling,  etc.  It  also  makes  it  possible  to  determine  the  draft, 
trim  and  heel  for  any  loading  condition. 

This  computer  program  was  developed  as  a tool  for  use  in  recreational 
boating  safety  research.  Certain  characteristics  of  recreational  boats 
made  a computer  program  of  this  type  desirable.  These  characteristics 
of  recreational  boats  are: 

a.  The  hull  shape  of  recreational  boats  (hard  chines,  cathedral 
hull  forms,  etc.)  precludes  usage  of  Simpson's  Rule  Integration  or  any 
other  parabolic  method. 

b.  The  low  length-to-beam  ratio  of  recreational  boats  causes  con- 
siderable coupling  between  restoring  forces.  A righting  arm  curve  which 
held  pitch  and  draft  constant  while  varying  the  heel  angle  would  have 
little  meaning  for  this  type  of  boat. 

c.  For  boating  safety  research  a program  was  needed  which  would 
give  equilibrium  characteristics  for  a given  load  condition  so  that  the 
effect  of  various  loading  configurations  could  be  examined. 

A general  consent  should  be  made  regarding  the  errors  Inherent  in 
representing  a hull  shape  by  a series  of  offsets.  If  too  few  offsets 
are  taken  to  define  the  hull  form,  errors  may  be  large;  while  If  too 
many  are  taken,  the  time  spent  in  preparing  the  data  Is  excessive.  The 
"trapezoidal  rule"  employed  for  integration  in  this  program  has  no  error 
by  itself — any  error  comes  from  Inadequacy  of  hull  representation.  The 
overall  error  depends  primarily  on  the  "engineering  judgment"  of  the 
user  in  preparation  of  the  offsets,  and  the  only  true  test  of  whether 
sufficient  offsets  have  been  used  is  comparison  with  full-scale  results 
for  the  particular  boat  in  question.  As  an  example  a sample  run  has 
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been  included  In  Appendix  C of  this  report.  For  the  equilibrium  case, 
results  were  within  the  accuracy  of  fullscale  measurements  of  actual 
draft,  trim  and  heel. 

The  program  and  subroutines  developed  are  written  in  FORTRAN  V language 
and  are  intended  for  batch  processing  use.  Section  2.0  of  this  report 
describes  the  various  sections  of  the  computer  program  and  the  subroutines. 
Section  3.0  provides  usage  instructions.  Section  4.0  explains  the  output 
from  the  program,  and  Section  5.0  explains  program  error  messages.  Included 
in  the  appendices  are  listings  and  flow  charts  of  the  program  elements, 
definitions  of  program  variables,  a sample  run,  and  test  results  used  for 
program  verification. 

2.0  PROGRAM  DESCRIPTION 

The  main  element  of  this  program  performs  the  input  and  output  and 
most  of  the  calculations  with  the  exception  of  the  integrations.  Figure 
1 shows  a block  diagram  of  the  program.  The  functions  of  the  various 
modules  shown  are  explained  in  the  following  sections. 

Three  subroutines  are  used  to  handle  the  integration.  Subroutine 
SECT  performs  integration  over  one  station;  subroutine  ENDPT  determines 
the  fore  and  aft  end  points  of  the  underbody;  and  subroutine  INTEG  integrates 
over  the  length  of  the  hull  (between  the  two  end  points).  These  subroutines 
are  discussed  in  Sections  2.10,  2.11  and  2.12  respectively. 

2.1  Section  for  Input  and  Output  of  Hull  Form  Data 

Input  and  output  device  numbers  are  defined  at  the  beginning 
of  the  program.  These  must  be  changed  for  each  computer  system.  This 
section  reads  in  data  which  defines  the  number  of  stations,  station  numbers, 
and  at  each  station  the  number  of  offsets,  the  longitudinal  position, 
and  the  half-breadth  and  height  at  each  point.  The  format  of  the  data 
cards  is  given  in  Section  3.3.  In  this  section  all  the  input  data  is 
also  printed  out  to  allow  a check  of  the  input  data. 

2.2  Section  for  Input  of  Case  Data 

In  this  section  the  number  of  cases  to  be  examined  is  read 
into  the  program.  The  program  then  begins  a DO  LOOP  in  which  it  reads 
in  data  for  one  case  and  then  does  the  calculations  for  that  case  before 
reading  in  the  data  for  the  next  case.  Within  this  section,  the  data 
for  one  case  is  read  in.  The  data  read  in  is  weight  of  the  craft,  LCG, 

TCG  and  VCG  (longitudinal,  transverse  and  vertical  centers  of  gravity) 
of  the  craft,  draft,  trim  angle,  heel  angle,  water  density  and  three 
control  codes  which  control  the  equilibrium  process.  The  format  of  these 
data  cards  is  given  in  Section  3.3.  Further  explanation  of  the  case 
data  is  provided  in  Section  3.2. 
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2.3  Beginning  of  Calculations 


Preliminary  calculations  are  performed  in  this  sc-tion.  Variables 
calculated  in  this  section  are  used  throughout  the  program.  This  section 
also  includes  the  DO  LOOP  which  controls  the  equilibrium  seeking  process. 

2.4  Transformation  of  Coordinates  to  Heeled  Coordinates 


This  section  rotates  the  center  of  gravity  coordinates  and 
the  offset  coordinates  from  the  input  coordinate  system  to  a heeled  coordi- 
nate system  which  is  perpendicular  to  the  water  surface  in  the  transverse 
plane  but  is  not  rotated  from  the  input  coordinates  in  the  longitudinal 
plane. 

2.5  Determination  of  Hydrostatics 

This  section  first  determines  the  water  level  at  each  station, 
which  is  a function  of  the  draft  and  the  trim  angle.  It  then  integrates 
the  individual  stations  by  calling  subroutine  SECT  once  for  each  station. 
This  determines  the  area  and  the  first  moments  of  the  area  in  two  direc- 
tions and  the  waterplane  beam  and  first  and  second  moments  of  the  waterplane 
beam.  Other  factors  determined  are  the  girth  of  the  station  to  the  water- 
line and  the  draft  at  that  particular  station.  After  the  station  charac- 
teristics are  all  calculated,  subroutine  ENDPT  is  called  to  determine 
where  the  forward  and  after  ends  of  the  waterplane  occur.  Next  subroutine 
INTEG  is  called  to  integrate  between  the  forward  and  after  ends  of  the 
waterplane.  Subroutine  INTEG  is  called  several  times  using  different 
parameters  as  the  argument.  From  this  is  determined  the  displaced  volume, 
the  longitudinal  center  of  buoyancy,  the  transverse  center  of  buoyancy, 
the  vertical  center  of  buoyancy,  the  waterplane  area,  the  longitudinal 
center  of  flotation,  the  longitudinal  waterplane  inertia,  the  transverse 
center  of  flotation,  the  transverse  waterplane  Inertia,  the  wetted  surface, 
the  metacentric  heights  and  the  righting  arms. 

2.6  Checking  for  Equilibrium 

In  this  section  a check  is  made  first  to  determine  whether  an 
unstable  equilibrium  condition  exists.  An  unstable  equilibrium  condition 
would  exist  if  the  waterplane  area  was  less  than  or  equal  to  zero  or 
the  Gil  (metacentric  height)  in  either  direction  was  less  than  or  equal 
to  zero.  If  one  of  these  conditions  Indicated  unstable  equilibrium  (i.e., 
negative  GM  or  no  waterplane)  in  a case  where  equilibrium  was  being  investi- 
gated, an  error  message  would  be  printed  out,  as  described  in  Section 
5.0.  Otherwise,  a calculation  of  the  residual  vertical  force  and  residual 
transverse  and  longitudinal  moments  is  made  to  determine  if  these  are 
within  predefined  error  limits.  If  the  three  parameters  are  outside 
the  error  limits,  then  new  calculations  of  the  draft,  trim  and  heel  are 
made,  and  the  sequence  is  repeated  for  those  new  conditions.  If,  however, 
the  three  parameters  are  within  the  error  limits,  equilibrium  has  been 
achieved  and  the  program  moves  on  to  the  supplementary  calculations. 

After  99  iterations,  if  equilibrium  has  not  yet  been  achieved,  a warning 
message  is  printed  out,  as  explained  in  Section  5.0. 
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2.7  Supplementary  Calculations 

This  section  determines  parameters  that  are  based  upon  those 
calculated  In  the  previous  section.  Parameters  which  are  calculated 
here  are  the  maximum  section  area,  the  maximum  draft,  the  maximum  beam, 
centers  of  buoyancy  and  flotation,  coefficients  relating  the  hull  volume 
and  vaterplane  area  to  overall  dimensions,  and  residual  vertical  forces 
and  righting  moments.  The  centers  are  then  rotated  back  to  the  original 
coordinate  system. 

2.8  Output  of  Results  by  Section 

This  section  prints  out  a table  of  results  by  section  in  the 
format  shown  in  Section  4.0.  The  information  included  in  this  table 
is  all  of  the  information  calculated  for  each  station.  At  the  end  of 
the  table,  the  case  data  which  was  used  is  printed  out  in  the  same  format 
in  which  it  was  read  from  the  data  card. 

2.9  Output  of  Complete  Results 

The  output  from  this  section  constitutes  the  second  page  of 
the  output  for  each  case.  This  is  a compact  table  which  presents  for 
each  case  the  input  data,  the  calculated  displacement,  centers  of  buoyancy, 
centers  of  flotation,  waterplane  area,  residual  moments  and  forces,  coefficients, 
stability  characteristics,  and  maximum  dimensions  of  the  underwater  hull 
shape.  This  output  is  further  described  in  Section  4.0. 

2.10  Subroutine  SECT  for  Calculating  Section  Characteristics 

This  subroutine  takes  the  offsets  for  one  station  and  integrates 
them  between  specified  limits.  First  and  second  moments  are  taken  using 
trapezoidal  rule.  Values  calculated  in  this  subroutine  include  section 
area,  section  beam,  section  waterplane  beam,  transverse  and  vertical 
moments  of  the  area,  transverse  first  and  second  moments  of  the  waterplane 
beam,  and  girth  of  the  underwater  sections.  , 

Subroutine  SECT  performs  the  integration  by  dividing  the  submerged 
section  up  into  vertical,  trapezoidal  shaped  strips,  each  defined  by 
two  consecutive  offsets.  The  area  of  each  strip  is  then  found,  along 
with  the  moment  of  that  area  about  the  centerline  and  waterline.  By 
sunning  up  all  the  strips,  the  area  and  center  of  the  station  are  determined. 

At  the  same  time  as  performing  this  integration,  the  subroutine  also 
determines  the  width  of  the  waterplane  at  the  station,  which  allows  the 
characteristics  of  the  waterplane  (area,  center  and  inertia)  to  be  determined 
as  well. 

2.11  Subroutine  ENDPT  for  Calculating  End  Points  of  Waterplane 

This  subroutine  takes  as  input  the  drafts  calculated  at  each 
station  and  determines  where  the  forward  and  after  ends  of  the  waterplane 
are  by  determining  where  the  drafts  change  from  positive  to  negative 
and  performing  a linear  interpolation  to  identify  the  crossing  point. 
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2.12  Subroutine  INTEG  for  Integrating  Over  Length  of  Hull 


This  subroutine  takes  as  input  X and  Y arrays,  where  the  X- 
array  is  the  longitudinal  location  and  the  Y-array  is  the  parameter  to 
be  integrated,  and  the  end  points  of  the  waterplane  defined  previously. 

The  subroutine  then  performs  a trapezoidal  integration  along  the  X axis 
between  the  two  end  points  to  determine  the  area  and  the  first  and  second 
moments  of  that  area  for  the  parameter  being  integrated. 

Subroutine  INTEG  performs  the  integration  by  dividing  the 
curve  (such  as  a "section  area  curve")  up  into  trapezoidal-shaped  strips, 
with  each  strip  extending  from  one  station  to  the  'ext.  The  area  of 
each  strip,  along  with  the  first  and  second  moment  of  that  area  about 
amidships,  is  then  found.  By  summing  all  the  strips,  the  area  under 
the  curve  (and  the  associated  first  and  second  moments  of  the  curve) 
is  determined. 

3.0  USAGE  INSTRUCTIONS 

The  following  sections  present  the  information  needed  to  run  the 
computer  program  in  its  many  different  forms. 

3.1  Hull  Form  Description 

The  axis  system  used  to  describe  the  hull  form  is  a right- 
handed  Cartesian  coordinate  system  with  the  center  of  the  system  located 
at  amidships  on  the  centerline  of  the  boat  at  the  keel.  X is  positive 
forward;  Y is  positive  to  port;  and  Z is  positive  upwards.  Angular  displace- 
ments follow  the  standard  sign  convention  used  on  right-handed  coordinates. 

The  hull  is  defined  by  an  arbitrary  number  of  sections  through 
the  hull  in  the  Y-Z  plane.  These  stations  may  be  located  at  arbitrary 
spacing  to  provide  the  best  hull  definition  in  areas  of  high  curvature. 

Figure  2 shows  station  locations  defining  two  typical  hu  1 forms.  Station 
numbers  or  names  are  also  arbitrary  and  bear  no  relation  ro  the  axis 
system.  A station  with  zero  area  may  be  used  at  the  bow  and  at  the  stern, 
if  desired,  to  allow  integration  to  continue  out  to  the  true  end  of  the 
boat.  The  program  is  currently  dimensioned  to  accept  up  to  25  stations. 

Each  station  is  defined  by  ar  arbitrary  number  of  points  located 
at  arbitrary  spacing  around  the  section.  The  station  is  assumed  to  be 
symmetrical  about  the  centerline.  This  assumption  of  symmetry  is  for 
convenience  in  input  of  hull  form  data  only — if  the  input  section  were 
changed,  non-symmetrical  hulls  could  be  accommodated.  Y and  Z coordinates 
define  the  shape  of  one  side  of  the  station,  and  the  program  creates 
offsets  for  the  other  side  of  the  station.  The  points  for  each  station 
are  entered  in  order,  as  shown  in  Sections  3.3  and  3.4,  with  the  first 
point  located  on  the  centerline  and  all  following  points  in  order  representing 
points  along  the  hull  bottom  and  up  the  sides  of  the  boat.  Figure  3 
shows  points  defining  several  typical  sections.  The  program  is  dimensioned 
to  allow  up  to  25  points  on  each  half  station.  (This  makes  50  points 
on  each  full  station  after  the  computer  creates  a mirror  image  about 
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FIGURE  2.  CHOOSING  STATION  LOCATIONS 
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the  centerline.)  The  only  restrictions  on  number  or  location  of  points 
are,  first,  that  at  least  one  point  must  be  given  for  each  station  (this 
would  define  a station  with  zero  area);  second,  the  first  point  must 
be  on  the  centerline;  and  third,  the  points  must  be  in  order  around  the 
hull. 

3.2  Condition  Information 

The  first  item  of  data  for  the  condition  information  must 
be  the  number  of  cases  to  be  considered.  This  can  be  any  integer  value 
greater  than  zero.  Each  of  the  cases  has  all  the  information  defining 
its  condition  and  the  type  of  analysis  to  be  done  on  one  data  card. 

The  information  included  on  this  card  is  weight  of  the  boat  in  pounds, 
longitudinal  transverse  and  vertical  centers  of  gravity  in  feet,  draft 
in  feet,  trim  and  heel  in  degrees,  water  density  in  pounds  per  cubic 
foot,  and  three  one-digit  control  codes  which  control  what  values  are 
to  be  allowed  to  reach  equilibrium.  The  different  methods  of  operation 
of  the  program  (which  are  controlled  by  the  control  codes)  make  use  of 
different  items  of  data  on  the  data  card. 

For  the  three  control  codes,  a zero  indicates  the  parameter 
in  question  is  free  to  seek  equilibrium,  while  a one  indicates  the  parameter 
is  held  constant.  The  first  digit  controls  the  draft,  the  second  digit 
controls  the  trim,  and  the  third  digit  controls  the  heel  angle. 

3.2.1  All  Parameters  Fixed 
(Control  Codes  111) 

In  this  mode  of  operation,  hydrostatic  properties  are  calculated 
for  one  condition  only.  That  condition  is  specified  by  the  draft,  trim 
and  heel  angle  that  are  put  in  on  the  condition  information  card.  Weight 
and  center  of  gravity  are  used  only  to  reference  metacentrlc  height  (GM) , 
righting  arm  (GZ) , and  restoring  forces  and  moments  to. 

3.2.2  Two  Parameters  Fixed,  One  Parameter  Free 
(Control  Codes  110,  101  or  Oil) 

These  three  modes  of  operation  are  used  to  investigate  restoring 
forces  and  moments.  For  Instance,  with  trim  angle  held  constant  at  zero 
degrees  and  draft  allowed  to  seek  equilibrium  conditions,  heel  angle 
could  be  varied  through  a range  of  angles  in  order  to  Investigate  the 
transverse  stability  of  the  boat.  In  this  way,  trim  angle  would  be  always 
equal  to  zero  degrees,  the  displacement  would  always  be  correct,  and 
the  heel  angle  would  be  the  specified  angle  in  the  different  cases. 

The  same  technique  can  be  used  to  obtain  pitch  righting  arm  curves  and 
heave  restoring  force  curves. 

3.2.3  One  Parameter  Fixed,  Two  Parameters  Free 
(Control  Codes  100,  010  and  001) 

This  mode  of  operation  is  useful  for  determining  the  restoring 
forces  and  moments  on  boats  like  recreational  boats.  For  example,  if 
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the  heel  angle  is  held  constant  at  each  of  several  angles  in  different 
cases,  while  both  pitch  and  draft  are  allowed  to  seek  equilibrium,  a 
transverse  righting  arm  curve  will  be  obtained  for  the  boat  in  which 
variations  in  trim  due  to  the  underwater  hull  shape  are  taken  into  account. 
This  can  also  be  used  to  determine  the  effect  of  heel  angle  on  trim  angle. 
In  these  modes  of  operation,  the  weight  of  the  boat  and  the  center  of 
gravity  of  the  boat  are  used  and  the  draft,  trim  and  heel  given  on  the 
condition  information  card  are  used  either  as  starting  conditions  or 
as  constant  conditions,  depending  on  which  parameters  are  allowed  to 
remain  free  and  which  are  fixed  conditions. 

3.2.4  All  Parameters  Free 
(Control  Codes  000) 

This  mode  of  operation  is  used  to  determine  the  equilibrium 
draft,  trim  and  heel  of  a given  hull  form  with  a specified  weight  and 
center  of  gravity.  In  this  case  the  draft,  trim  and  heel  values  given 
on  the  condition  information  card  are  used  as  starting  values  for  the 
iteration. 


3.3  Data  Card  Format 

Hull  form  data  and  condition  information  data  are  entered 
on  standard  80-column  computer  punch  cards  in  the  format  indicated  in 
this  section. 


3.3.1  Hull  Form  Data  Format 

Card  number  one  of  the  hull  form  data  is  an  alphanumeric  title 
describing  the  boat  in  the  first  69  characters  on  the  card.  Card  number 
two  contains  the  number  of  stations  as  an  Integer  right-justified  in 
columns  1 through  10.  Card  number  three  contains  alphanumeric  station 
titles  or  numbers  in  a format  of  IX,  15A5.  If  more  than  15  stations 
are  used,  additional  cards  are  used  for  the  remaining  stations.  For 
the  first  station  from  the  bow,  a card  is  given  which  has  the  number 
of  points  on  the  station  and  the  X-distance  of  the  station.  The  format 
of  this  card  is  110,  F10.4.  Following  this  card  is  one  card  for  each 
point  on  the  station.  Each  of  these  offset  cards  has  two  numbers  on 
it,  the  Y-coordinate  and  the  Z-coordinate  for  the  point.  The  format 
for  these  cards  if  2F10.4.  This  sequence  is  repeated  for  each  station 
on  the  boat. 


3.3.2  Condition  Information  Card  Format 


The  first  card  of  this  section  contains  the  number  of  cases 
as  a right-justified  Integer  in  columns  1 through  10.  The  remainder 
of  that  card  can  be  used  for  description,  if  desired,  but  this  description 
is  not  read  by  the  program.  Following  the  first  card  in  this  section 
are  a number  of  condition  information  cards  equal  to  the  number  of  cases. 
Each  oi  these  cards  contains  the  information  of  weight  in  pounds,  longitudi- 
nal center  of  gravity  in  feet,  transverse  center  of  gravity  in  feet, 
vertical  center  of  gravity  in  feet,  draft  in  feet,  trim  and  heel  in  degrees, 
rho  (the  weight  density  of  the  water)  in  pounds  per  cubic  foot,  and  the 
three  control  codes.  These  parameters  are  formatted  on  the  card  in  a 
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format  of  IX,  8F9.3,  3X,  312.  Figure  4 shows  a condition  information 
card  coding  form,  which  is  helpful  in  preparing  the  data  for  keypunching. 

3.4  Sample  Input 

Figure  5 shows  a sample  input  for  the  hull  form  of  a jon  boat. 
Figure  6 shows  the  condition  information  for  this  same  jon  boat. 

4.0  INTERPRETING  THE  OUTPUT 

The  first  section  of  the  output  is  a repetition  of  the  hull  form 
data  used.  The  data  which  defines  the  boat  hull  form  is  printed  out 
in  the  same  format  as  the  cards  with  the  input  values.  This  allows  checking 
of  the  hull  form  data.  The  total  number  of  cases  is  also  printed  out 
in  this  section.  See  Appendix  C for  an  example. 

The  results  for  each  case  appear  on  two  pages  of  output.  The  first 
page  of  the  output  contains  the  results  of  the  intermediate  calculations 
for  each  station.  The  first  line  of  the  output  gives  the  draft,  trim 
and  heel  in  original  coordinates.  Directly  under  the  draft  and  trim 
appear  the  values  in  the  rotated  coordinate  system.  Appendix  C shows 
a sample  of  this  page  of  output.  Directly  under  this  table  is  presented 
the  case  data  used.  This  is  the  same  data  that  appears  on  the  condition 
information  card  in  the  same  format. 

The  second  page  of  output  for  each  case  is  a presentation  of  the 
overall  calculations.  Appendix  C presents  a sample  of  this  output. 

In  this  output  values  that  appear  in  parenthesis  are  in  the  original 

coordinate  system,  while  values  not  in  parenthesis  are  in  the  rotated 

coordinate  system.  Where  only  one  value  is  given,  both  coordinate  systems 

correspond.  Block  number  one  of  this  output  presents  the  draft,  trim 

and  heel  of  the  boat  as  calculated,  as  well  as  the  weight  and  center 

of  gravity.  Block  number  two  presents  the  calculated  displacement,  center 

of  buoyancy,  waterplane  area,  and  center  of  flotation.  Block  number 

three  presents  overall  waterplane  length,  waterplane  beam,  maximum  draft, 

maximum  section  area,  the  location  of  the  ends  of  the  waterplane,  and 

the  wetted  surface.  Block  number  four  presents  the  residual  moments 

and  forces.  (In  an  equilibrium  case,  these  forces  and  moments  would 

be  zero,  or  very  nearly  so.  Equilibrium  is  considered  to  be  reached 

when  changes  of  less  than  0.001  ft.  in  draft  and  0.0001  radians  in  the 

two  angles  are  called  for.)  Block  number  five  presents  stability  information, 

both  longitudinal  and  transverse;  and  block  number  six  presents  coefficients. 

5.0  PROGRAM  ERROR  MESSAGES 

The  program  includes  one  warning  message  and  one  error  message. 

The  warning  message  appears  in  any  case  where  equilibrium  is  being  sought 
and  no  equilibrium  has  been  reached  after  99  iterations.  The  format 
of  this  warning  statement  will  be: 

"WARNING-NO  EQUILIBRIUM  XX  XX  .XXXXXEtXX  .XXXXXEtXX  .XXXXXEtXX" 
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FIGURE  5.  SAMPLE  OF  HULL  FORM  INPUT 
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FIGURE  6.  SAMPLE  OF  CONOITION  INFORMATION  INPUT 


The  first  number  after  the  warning  indicates  the  number  of  Iterations 
attempted.  The  second  number  indicates  the  total  number  of  equilibrium 
conditions  which  were  not  met.  This  can  be  either  one,  two  or  three. 

The  last  three  numbers  in  scientific  notation  form  are  the  change  in 
draft,  change  in  trim  and  change  in  heel  angle  from  the  last  Iteration. 

The  second  diagnostic  generated  by  the  program  is  of  the  form: 

"CALCS. STOPPED-UNSTABLE  .XXXXXE±XX XXXXXE+XX  .XXXXXE+XX" 

The  numbers  following  this  error  message  are  draft,  trim,  heel, 
waterplane  area,  longitudinal  metacentric  height  (GML)  and  transverse 
metacentric  height  (GM_) . This  error  and  message  occurs  when  one  of 
the  restoring  forces  is  too  small  for  stability. 

While  neither  of  these  error  messages  prevents  further  execution 
of  the  program,  their  appearance  in  the  output  causes  the  results  to 
be  suspect  for  that  particular  case. 
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APPENDIX  A - PROGRAM  LISTINGS  AND  PLOW  CHARTS 


Lilting  of  Main  Program 


HYDROSTATICS  and  equilibrium  characteristics  program 

DEVELOPED  FOR  RECREATIONAL  BOATING  SAFETY  PROJECT 
DEVELOPED  AT  USCG  R*D  CENTER  GPOTONzCONN. 

OEVELCPEO  EY  PETER  SECFEST  SEPTEMBER# 1975 

INTEGER  TITLEC23)/NPTSC25) 

REAL  XC25>/YC25z5Q)/ZC25z50>/Y0(25/5U>zZ0C25z50)zNUMaCRC2S  ) 

REAL  LCGzLCG0zZL0MC25)zAC25)/AYMC25)zAZMC25)zfcC25)zbYMC25> 

REAL  JYI(25)/3MPC2S)/gIRTHC2S)/CRaFTc2S)/LCBM/LC3/LCFH/LCF 
REAL  LMP I/CA1 C25  J/CA2C25  >z*YCC25  >z AZCC25>zBYCC25>zLPCMzLML 
C SECTION  FOR  INPUT  AND  OUTPUT  OF  HULL  FORM  DATA 
INPUT®? 

IPR I NT=4 

MR ITcC JPRINT/42) 

NEADCINFUTzUCT ITLE(J)zJ=1z23) 

MR ITECIPP INT/1 ) (TITLE! J)z J=1 *23  ) 

READ  < INPUT  / 2 )NS  T A 

MR  ITECIPRIMzfcJKSTA 

READ(  INPUT  #5  XNUM3EKJ  ) z J*  1 z N 3 T A ) 

MR  I TEC IFR INTz5)C NUMBER C J >zj*1 zNSTA  ) 

00  100  N1=1#NSTA 

READC INPUT/2)  NO  FFSzXCM  ) 

MR  I TEC IPR IF Tzo)NOFFSzX CN1 J 

NPTS  CN1 )®2»R0FFS 

00  100  R2*1 /NUFFS 

REAOC INPUT  z4  )YOC  N1  zN2  > zZOCM  zN2) 

MR  ITECIPR  IRT/7)Y0CN1  / h 2 )/Z0CN'1/N2) 

NJ*NPTSCN1  ) -N  2*  1 
Y0CN1/N2)®-Y0CN1 zN2) 

Z0CN1/N3 )*Z0CN1/N2) 

IuO  CONTINUE 

C btCI ION  FOP  INPUT  OF  CASE  DATA 
FEAD  C INPUT  z t INCASES 
MR ITtCIPRINT/UTNCASES 
CO  2uC  M»1/NCASES 

READCINFJTzSIMEI GHT/LCGUzTCCO/VCGO /Tu/THLT AUz PHIO/PF 0/ NT/NT HE TAz 
iNPHI 

C BEGINNING  CF  CALCULATIONS 
VOLU  ME  = ME  IGHT/PHO 
PhI*FHIO/S7.2957S 
CUSPHI-COSCPHI) 

S1NPHI-SINCFHI) 

T»TO»COSPHI 

TFiETa  = THETAL»C0SPHI/57»2957E 

LCG*LCGO 

00  3U0  I *1 z9V 

C TRANSFORMATION  OF  COORDINATES  TO  HEELED  COORDINATES 
TCG=TCGO»CCSPHI -VCGO»SI NPHI 
VCG® VCGOaCCSPHIzTCGOaSI NPHI 
00  350  1 1 ®1 zNSTA 
ZLOM  C 11 )**1G0« 


A-l 


I P TS  =NP  T S C 1 1 ) 

DO  35J  I 2*1 . IPTS 

V<I1.I2)*YG(I1»I2>*COSPhl-ZO(I1»I2)«SlKPHI 
Z<  II. I2)=/L(I1.I2)»COSPHl*YO(I1.I2  )»S1NPHI 
IF<ZUl.I2).LT.ZLUV(I1))ZCOIKI1)*2<I1,I2> 
iaO  CONTINUE 

C DETl  kMNATION  OF  HYCROSTATICS 
00  4 UO  K 1 =1 . NS T A 
ZHIGH*T*X(K1)«TAN(THETA) 

NP  = NPTS(K1  ) 

CALL  SECT  (K1.NP.Y.Z.ZH1GH.A1  . A Y M 1 . A 7M 1 . fi 1 . 6 Y M 1 » BY  1 1 » B UP  1 . u I r T Hi  ) 

A(K1  ) = A 1 

A Y H ( K 1 ) = A Y **1 

AZM(K1)=AZP1 

B ( K 1 ) =B  1 

BYNUM  > = BYH1 

BY  I ( K1 )*BY 11 

bUPCKl )*bUP1 

G I R T H (K  1 ) *6  I F 1 HI 

DRAFT (KI)-ZHIGH-ZLOU(Kl) 

4oC  CONTINUE 

CALL  LNDPKI.NSTA.X. DRAFT. XSTAAT.ISTaRT) 

CALL  cNDPTOSTA.NSTA.X.OPAFT.XEND.ILNU) 

IF < I END.LT. ISTAK T.CF .NSTA.LE .1 )U0  TO  410 

CALL  INTECU.A.XSTART.  XE  ND»  ISTAP  T.  IENO.VDISF  .LCBH.EL  ANK  ) 

CALL  INTt&(X.AYH.XSTART.XEN0.1START. It NC.TCBM.bLANK. BLANK) 

CALL  INTEGt X.AZH.XS TAR T.XE NO. IS  TAR  1 .IcND.KCeM. BLANK .BLANK  ) 

CALL  INTEG<X»3.XSTAFTiXENQ.lSTAKT»IEND»AuP»LLFh*LUFl) 

CALL  IN TEGtX. BY*. XSTART.XENC.lSTART.lt NO. TCFM. CLANK. BLANK) 

CALL  INTEGCX. iY I. XSTART.XcNC.ISTART.il  NO. TUF I. BLANK. BLANK) 

CALL  INTE»(X. GIRTH. XSTART.XEND. 1ST *CT»IFND»U3.PLANK, BLANK) 

GO  TO  42 J 
410  CONTINUE 
VD ISP=0. 

LCBM*u. 

TC6**  = J. 

VCBM*P. 

A U P = u • 

LCFR*C. 

LUP 1*0. 

TCFF*u. 

TUP  l *0. 

US*G  • 

4*0  CONTINUE 
LCF*0. 

T CF  *«• 

I F ( A UP.LL .0 .0001  )G0  TO  4.JU 
LCF*LCFH/AGP 

tcf*icfp/ave 

L W P I *LWP  I-l CFM.L CF 

TBPI-TBF i-tcfnrtcf 

4 iC  CONTINUE 
LCB-C. 

T CB*u  . 

V CB*Q  . 

BML*U . 

8PT*u. 


A-2 


W 


IF (VDISP. Lb. O.OOOI  )G0  TO  410 

LCBmLCBH/VD  ISP 

TCe-TCBH/VDISP 

VC6-VCBM/VG ISP 

BHL*LUP I/VO  ISP 

BNT  = TUP  I / YD  I S? 

440  CONTINUE 

GHLxbHL*YCB-YCG 

fcPTxbMT»vce-vct 

G2L=(lCC-LCfc)«C0S(THETA)»(VCG-VCB>»SIN(THtTA) 
GZT=TCB-TCL 
LWLxXSTART-XENO 
C CHECKING  FOR  tQUILIbRIUP 

IF (NT  .LE.0.AND.AUP.LE.C.0G1  )G0  TO  500 
IF (NTHETA.LE.O.AND.GNL.Lt.O.OUl  )G0  TO  SOO 
IF (NPHI .LE.0.AND.GHT.LE.U.001  )G0  TO  500 
I TE  S T»0 
CT«0. 

OTHETA-O. 

OPHI >0* 

IF (NT.GT .0)40  TO  510 
0Tx( VOLUHE-VOISF J/AWP 
IF(AbS(0T).GT.0.001)ITEST*=ITE£T*1 
510  CONTINUE 

IF(NTHETA.CT.O)GO  TO  520 
DTHE  T A*G  ZL /UHL 

IF(AbS(DTHETA).CT.0.0uC1)ITEST=ITEST»1 
5c0  CONTINUE 

IF ( NPH I • GT *0  ) GO  TO  530 
OPHI«GZT/GHT 

I F ( A 8S( DPH I ) »LT .0.CGU1  )I TEST ■ I TE  ST ♦ 1 
530  CONTINUE 

IF ( I T tS  T »LE  *0  )G0  TO  540 
T*T*OT-TCF«DPHI*LCF "OThETA 
ThET  A »(  THcTA  + DTHcTA  ) / C OC Prl  I 
P H I XP  HI ♦ OPH I 
COSPHIxCOS(PHI) 

S I NPh  I*S  I N ( PH  1 ) 

THETA *THETA»COSFHI 
IF (NT.GT. 0)T«T0»C0SPrtl 

IF  (NTHF  T A.GT .Q)THETAxT  HET  AOaCOSPhI /57.2G57E 
w oO  CONT INUE 

NR  ITEdPMNT.IOTI.  I TfcS  T »Ul  .0  T HE  TA  . DPH  I 
GO  TO  54  J 
500  CONTINUE 

GRlTt(IPRINTf11)T.THETAfPHWAUP«GHL>GHT 
54C  CONTINUE 

C GuFPLE  ME  N IAR  Y CALCULATIONS 
A h A X xu. 

OR  N I X s0 . 

BEAR-U. 

IF  (IENC.LT.ISTART.CF.NSTA.lt  .1  )G0  TO  o2C 
00  605  J-1.NSTA 
ATC(  J)-0. 

AZC(  J)-0. 

B Y C ( J )-0  . 

CA1 ( j)»0. 


CA  2 ( J)»C. 

6u5  CONTINUE 

00  60U  J<ISTAI<T<IENC 
IF  (A(J)»GT«AMAX  )AHAX«A(J) 

IF(D*AFT(J>«CT«DMAX>DPHAX*0RAFT(J> 

IF(BWPCJ>.GT.8EAM)6EAH«eBP<J> 

6uC  CONTINUE 

DC  6 1 L J< 1ST  ANT , ItNC 
1F(A(J).LE.C.J0C1 )G0  TO  6lC 
ATC( J)=AYN(J)/A( J) 

AiCt J)<AZH(J)/A< J) 

IF  <fc< J).GT.C.0001  )P1C(J)<6TN(J)/BCJ) 

CAl <JJ*A(J)/AMAX 

IF <9WP( J ).LC .JtOOOl .OR .DRAFT ( J) .LE »u«UP01 )GC  TO  610 
CA2(-)<»CJ)/BWFCJ)/CRAFT(J) 

610  CONTINUE 
e>2C  CONTINUE 

ThETaD=THET>»57. 2957a 

PHIO*PHI »57. 29578 

D I SP  <V0 I SP  *A  rtU 

LMON*JISP»C.ZL 

TN0H<0ISP«G7T 

VFOr  tt*C  ISF  -BE  IC-HT 

IF  (L  BL»LE  »l  .U  J1  .OF  . tE  A M .LL  »0  .001  .OF  .DAHAX  «l  t .0  .001  ) C-  0 TO  o3U 

Cb<VuISP/LWL/BEAl/rF.iAX 

CX=AHAX/BEXN/,DAHAX 

CMCB/CX 

CH=A  BP/LWL/EtAN 

C VP  = CB/C  B 

CIL=LnPI/LGL*»3/aEAN 
CIT  = TWPI/LWl/6EAi1»»3 
GO  TO  64 J 
o JO  CONTINUE 
CB  = 0 • 

C X < U • 

CP*U. 

CB=L  • 

CVP<U. 

CIL=U. 

Cl T=u. 

c 4C  CONTINUE 

T1=T/C0SPHI 
ThEl A1<THETAU/C0SPHI 
TCG1=TCG»CCSPHI*VCG«S1NPH1 
VCG1*VCG«C0SPHI-TCG»SINPHI 
TCS1*TCE«C0SPHl*VCP»SlNPMl 
VCB1 <VCB»COSPmI -T  CB»S INPHl 
TCF1 * T CF  »COSPHI*T  »S 1NPHI 
VCF1=T«C0SPNI-TCF«SINPHI 
C OUTPUT  OF  AESULTS  BY  SECTION 
BM  TE  ( IPR  INT<40  ) 

bP1TE(IP'<INT<52)<TITL£<J)<J=1»23) 

bP  ITfc  (IPRINT»20)(T1 # THE  TAl < P H 1 D»  T < THE  T A D ) 

NS=N3T*/10» 

NS1 =1 U»NS 

IF (NSI.l T.NSTA)NS=NS*1 
00  100  I S T A « 1 »NS 


J 
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L = I S T A» 1 3 
K ®L  - V 

IF (L.GT.NSTa >L=NSTA 
WRITEdPRINT#21  )( NUMBER (J)»J=K»L  ) 

WRITEdFRIRT»22)(X(J).»JxK*L) 

WRITEdPRINTj23>(A(J)»JxK»L> 

WR  1TE  dPR  IR  T#24  > (CAU  #L  > 

WRITEdPRINT#25)(CA2CJ)#JxK,L> 

UPITE(IPRINT«26)(AYN(J)«J>K#L) 

WRlTtdPRlM»27XAYt(J>#J=KjL) 

WR  I Tl(  I PR  INT»2E  ) ( AZML  J )#J  = K *L  ) 

WRITEdPRlNT»29UAZCLJ)»JxK,L> 

WRITLdPRlN'Tj30)(BCJ)*JxK#L) 

MRITC(IFRIMjJ1)(8YP(J)*J  = K#L) 

WRITtdPKINT*  32  )(BYCCJ)/J=K#L) 

WR  ITE(IFRIKT»33)(BY  I(J)»J=K#L  ) 

WP  ITEdPRlNT#34MBWPt  > 

MRITE<IPRirT*35XDRAFT(J)»J  = K#L) 

WRITEdPRINT*36)(GIRTH<JXJ  = K,L) 

WR  I T t ( IPR  IM  *41  > 

800  CONTINUE 

WR  I T E d P R I N T » 4 5 ) 

WR  ITECIPRINT»9)HEICHT»LC(.0f  TCGO^VCtfO*  TO#THETAJ»PMlCfhiiCf  NT  .VTriFT  Aj 
JNPH  I 

C OUTPUT  OF  COMPLETE  RESULTS 
WP  I TEUPR  1NTj40> 

WRITC<IPMNT»aO) 

WR ITE(IPRINT*51  ) 

WRITt(IPMNT»52)(TITLE(  J)*J*1#23) 

WRITEdPRIRT»S1  ) 

WRITEdPRINT»50) 

WR  ITEdPRIM»S4)T1iTjl<ISP 

WR ITl(IPRINT»S5)ThETA1 »ThtTAD»VDISP 

WRITE(IPRINTj36)PMIC»LCB 

W R I T t ( I P R INT#S7)WEIGhT » T C 3 1 * T C b 

WRITt(IPRINT»b8)LCGiVCb1#VCE 

WR  ITtdPRIR  T»:>9)TCC1  jTCGiAWP 

WR  I TEC  IPR  IN'T#60  ) VCG1  * V CC-»LCF 

WR  ITE(IPRIMT«61 JTCF1  »TCF 

WR  I T fc  < I P R If«T»62)VCF1  # T 

WRITc(IPRIKT#S3) 

WR  ITEdPRINT*  50) 

WR  ITEdPR  IRT»63)LWL»LM0M 
WRITE(IPRINT»64)BEAM#GZL 
WR  ITfcdPR  IPT*65)0RhAX 
WR  I T E ( I PR  IE:T  »6fc)AMAXfTM0M 
WRITELIPRIKT,67)XSTART*GZT 
WR  ITEUPRIR!T#68)XEN0 
WRITEd  PRINT*  69  )WS*V  FORCE 
WR ITtdPKIR'T*53) 

WRITE(IPRINT»S3) 

wRITE(IPRINT»SO) 

WR  ITEdPR  INT*  70)LWF  I*C6 
WR I TL ( IPR IN  T *71 ) 8ML  * CX 
WR ITEt IPR INT*72 >SNL*CP 
WR  ITEdPRINT*73)CW 
WR  lTtdPRINT»74)TWP  l » CVP 
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UR  lTfc(IPRlNT.75)BMT.ClL 
UR  ITEtIPRiNT»76)GHT.ClT 
UR  I T E ( IPR I N T » 53  ) 

UR  ITE(IFRINT»52) 

UR ITEUPRINT.50) 

UR ITEIIPRINT.77) 

UR  I Tt (IPR INT.50  ) 

2GC  CONTINUE 

1 FOPMAU23A3) 

2 FORMAT<5I1U) 

3 FORMA  T( 1 1Q.F10.4  ) 

4 F0RHAKSF1C.4  ) 

5 FORMATd  X.1  5A5  ) 

e FGRMATdX.IIG.F 10.4) 

7 FORMATd  X.2F1U.4  ) 

8 FCRHAK17H  NUMBER  OF  CASES=»I2) 

9 F0F.MATdX»fF9.3»1X.2I2  ) 

10  F0RMAT(23H  W ANN  I NG-  N 0 E Q U I L I BR  I UM  , 2 1 4 . 3E  1 2 • 5 ) 

11  FORMA  T( 23H  CALCS.ST0PPED-UNSTAELE.6E12.5  ) 

20  F0NMAT(///19H  RESULTS  BY  SE C T 10 N . 9 X ,c HT  = > F 7 • 3 . 1 2 H FT.  THE  T A = » 
SF7.3.13H  OEt.  PH1=.F7.3»5H  DEG..4X.2GH»»VALUES  IF  ORIGINAL. 

s 1 4 h CCORO I NATES*  »./20X  »F7.3. 1 3X. F 7 .3. 16X .4H.0JP.9X .1 1 H»  «VALUtS  IN. 
123F.  NOTATED  COORDINATES"*) 

21  F0RMAK1  X.2CHST  AT  ION  NUMBER  »10(3X.A5)) 

22  FORMATd  X.20HX-P0S  IT  ION  l F T ) , 1 OF  1 0 . ? ) 

22  FORMA T< 1 X .2CHARE A (SO .FT ) . 1 CF 1C .3  ) 

24  F ORM AT( 1 X .20HCA 1 =ARE A/MA X .AN E A .10F10.3) 

25  F0RMAT(1X»2CHCA2=AREA/8EAM/DRAFT  .1CF10.3) 
t6  F0RMAT(1X.2CHAREA  Y -MO ME  NT ( CU .F T ) . 1 0F  1 0 . 3 ) 

27  F0RMA7<  1 X »2ChAREA  Y-CENTEN  ( F T ) . 1 r F 1 C . 2 ) 

11  FORMATd  X.2CHAREA  7 - MO  ME  N T ( CU  *E  T ) . 1 Of  1 0 . 3 ) 

4.9  FORMATd  X.20HAREA  Z-CENTER  ( F T ) . 1 CF  1 0 . 3 ) 

30  FORMAT (1 X.20HUP  UIDTH  ( F T ) . 1 CF 10 .2  ) 

31  FORMAT (1 X.20HUP  Y-MOMENT  ( S 0 .F T ) . 1 CF 1 0 .3  ) 

32  FORMATd  X.2CHUF  Y-CENTER  (FT)#10F10.3) 

3?  FORMATdX.ZCHUP  Y-INERT  IA  ( CU .F T ) . 1 OF  1 0 . 2 ) 

34  F0R“AT(1X»4',CHUF  BEAM(OVERALL  )(FT  ).1CF1C»3) 

25  FORMATd X.2LHDRAFT  ( F T ) . 1 0F  1 0 . 2 ) 

36  FORMATd  X.20HGIRTH  ( F T ) . 1 CF  1 0 . 3 ) 

4C  FORMATd  HI  ) 

41  FORMAT(lHO) 

4.2  F0RMAT(2lHlHULL  FORM  DATA  USED-) 

42  FORMATd IHOIATA  USED-) 

50  F0RMAT(lX»72(1h»)) 

al  FORMA T( 1 X .1 R ». 70X.1 h"  ) 

52  FORMATd X .1 h".23A3»2H  ») 

53  F0RMaT(1X.1H".34X»1H"»35X»1H») 

54  FORMATCIX# 8H«T  ( . F 1 0 • 3 . 1 h ) , F 1 0 .5 . 2 5 H FT  "DISPLACEMENT  . 

JF 1 0. 3»8h  LE  «) 

55  FORMATd  X.bH»THE  TA  ( . F 1 0 • 3 » 1 H ) . F 1 0 . 3. 2 5 H DEC  "VOLUME  DISP.  » 

IF  10.3. 8H  CU.FT  *) 


56 

FORMATd  X.fcH«FFiI 
IBM  FT  «) 

(.F10.3.lH)»5X.1  SHU. 000 

Dt  G 

•LCE.15X.F10. 3. 

57 

F0RmaT(1X.19H"UE  IGHT  .F10.2.14H 

LB 

"TCE  ( »F  1U.3. 

i 1 H ) . F 1 0 . 3 . fe  H FT 

« ) 

58 

F0RMATdX.19M"LCG 

.F10.3.14H 

FT 

• VCb  ( .F  10.3. 

JlH)»F10.3.Zh  ft 

« ) 
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59  FOPMATdX.SF.ATCG 
1F1U.3.8H  SO.FT  ») 

(»F10.2.1H)»F10,!.25H 

FT 

•UATEFPLANE  AREA 

$ 

oC  FCPMATd  X.FHAVCC 
SF1C.3.8H  FT  a) 

(»F1G.3.1H).F10.3.25H 

FT 

• LCf 

f 

tl  F0PNATdX»lPA.34A»fcHATCF  C »F  1 C . 3 » 1 H ) . F 1 0 • 3. 8 H FT  • ) 

62  FORK  A T ( 1X»1h»»34X»SHAVCF  ( »F  1 0 • 3 » 1 h ) »F  1 o .3  » SH  FT  •) 

63  F0PNATdX»19HAWATEPPLANE  LfcNGTH  .FIO.S.I&H  FT  a L C N l»  .M  OML  U T » 7 X , 

1F1G.4.8H  FT. LB  a) 

64  F O FT  A*  A T £ 2i)H  awATEPPLANE  BE  A M .F10.3.25H  FT  «LONt.G2 
XF1G.5.8H  FT  ■ ) 

65  FGPNATd  X . 1 9rt»M A X I HUN  CRAFT  .F13.3a7h  FT  »»35X»1H«) 

66  F 0 0 N A T ( 2 JH  » M A X .SECT  ION  A • E A .F10»3#25H  SO. FT  aTRANS..1C*ENT 
XF1U. J.8H  FT .15  a ) 


67  FOFPAT<2‘jH  •FWO  WATEPPLANE  E NO  * F 1 3 . 2 5b  FT  aTRANS.„?  . 

!F10.3.Eh  FT  • ) 

tfc  FOPHAT(2JH  .AFT  WATFRPl ANE  E NC » F 1 3 . 3 . 7M  FT  a.JsX.Iha) 

69  F OR  P A T( 23H  » a E T T t G SURFACE  aF10.3»25H  SQ.FT  »\EF  T .CAl  FOrtc  . 
4F1G.3.8H  Lt  ■ > 

7C  FGrMAT(2GH  ALONG.  INEPT  IA  .F10.3.25H  FT4  »FLOCK  COEFFICIENT  , 

A4X#F3.5»5Xj1ri«) 

71  FOKHAT(1Dm  «B«L  *1GX #F1C »J»2fcH  FT  «*A X.SE CT  IGF  CGEF  . 

ilx .Fa. S.SX.Iha) 

72  FOPNATdjH  «GHL  » 1 CX  »r  10  .3  .2  6H  FT  «PA.I31AT1C  COLF  . 

43X.F3»5.5A.1nA) 

73  F0P**AT(1X#1f«»34X#2CH»WATEAPLANE  COEF.  . 3 X . F 2 .5  * 5 X » 1 h»  ) 

74  F0»PAT(1srt  ATPANS.INthTlA.3X.  FIG.  3»?6n  F T 4 » V E f T .F  F 1 S*  A T i C CGEF. 

4 j X .F  a.5.5X.1ri«  ) 


75  FOFNATdOH  »&1T  .10X.F1L.3.26H  FT  aL  C NG . I NE  N T I A COEF  > 

A 3 X . F a.S.SX.lHA) 

76  FOPNATdOH  »&MT  .1CX.F1 0.3.26^  FT  aTPANS  • IN'ChT  IA  COEF  . 

A 3X.F a.S.SX.lHA) 

77  F0FNATdX.1hA.4lH  VALUES  IN  PARENTHESIS  AhE  IN  NON-F  OTATfcu  . 

X 2 1 H ORIGINAL  C0CR01NATES.3X.lHA) 

STOP 

tNl> 
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Flowchart  of  Main  Program 
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I T»»»bfOHI»I  ION  OF  COCFOlMTkS 
• TO  HtElEC.  COOFDIIUTES 


' 


\ 
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a 


CALL  !HlEG(«»*t*»»ST4l  T ,l£  / 1SI  » ► I . 

imO»VCBf»^L»K»lL**>»  > 


CO  * T I WUE 
ICF-J. 
TCF«u. 


A-H 


G«T«BMT*VCB-¥CG 

W2l*(LCG-LCB  )«COS( THETA >«<VCG-VCE)»SIM 
T ME  T A ) 

C2T-TCB-TCG 

LML*ISTH»T-XEFP 


CHfcCK  IKU  FOR  EOU11IBF  IOH 


IMJI 


if  utmu.ii.cK 
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1*T*CT  -TCr  •OFHlUCr«OTHE  TA 
T*fcrA«<fH£  ta*o theta ) /cos phi 
Fh l«PM I«OFHI 


CO  NT 1NUE 


M«m(lPR  INT«11  )T«Th£TA#Ph1 
CHT 


IF(A ( J ).&T  .AHAI  ) 


{ r >¥«I VMV 


JD  JO  OH>(Ol8(DS(Bll)(l)ll)ll)IXIl8 


IF(B(J)«GT*G«0301  ) 
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GC  TO  640 


I 


A-21 


V 


* $»NS  JA/10 
N$1 *1 J«KS 


Ml  SI 
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TE  ( IPM 


»»lTMlPMNT»S5)THeT»1»THET*D»*0ISP 


r 


a 
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MR  ITE  < IPP  |M  lb  3)IML  »l  «0H 
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V*ITE<lPP[kr«72>tttl»CP 


**I  Tfc  < !P«  IM»73>Cm 


A-27 


it  F-«k*Tc23m  allXaSC CT !0»  Aft* 
F13a.4«2*4  $ ' J • F T •T*#N$#HC*t  NT  , 


I 


I 


i I I I I I I I I I I I I I 


I I I I I I I I I I I I I I I I I I 
• aw  » x « 

• t>  I • K a-O  X «»  • 

unu.  r-  % cm  *o«a  ‘-•ox 


>» 


h i 


Listing  of  Subroutine  SECT 


SUfc*QUTINE  SECT ( I . N P # Y tl i ZHI G H. A . A YM.AZM.EiBYM  .BY I . E W P . G IR  TH) 
DIMENSION  r(25.50)»Z(25.50) 

A = C • 

A Y N=  O. 

A Z N = U • 

0 = 0 • 
b YM=0. 

BY  1 = 0. 

GIPTH=0. 
f HAX--1G0. 

YNIN=t100. 

VI *Y ( IjI  ) 

Z1=ZHIGH-Z(  1.1  ) 

00  10  J=2fNF 
Y2  = Y ( I. J) 

D Y = Y c -Y  1 
Z2=ZH1GH-Z( I . J) 

D Z = Z 2-Z 1 

IF < Z 1 .GT  .0.0001  .AN0.Z2.GT. 0.0001  )G0  TO  11 
IF <Z1  .GT.O.C001 . AND. Z2.LE.0. CQ01  >G0  TO  12 
IF (Z1.LE.0.C001.ANC.Z2.GT.0.0001  )G0  TO  13 
GC  TO  14 

11  CONTINUE 

0 A =D  Y « ( Z 1 ♦ Z 2 ) /2  • 

YA VG=(Y1 ♦ Y2  )/2. 

DAY;‘  = Z1«OY»YAVG»DZ»DYMY1*2.«DY/3.>/2e 
DAZM=Z1»DY»Z1/2.*DZ«DY»(Z1tOZ/3.)/2. 

DGIRTh*SQRT (DY»CY*GZ»DZ) 

OC  TO  15 

12  CONTINUE 
DY=-DY«Z1/DZ 
0A=DY«Z1 /2» 

YA  Vt  * Y1 ♦OY/2 • 

DAYN»DA«<Y1*0Y/3.) 

D A ZM  *DA  « Z1 / 3 • 

Y 0 = V 1 *DY 

IFC YO.GT.VMAX)YMAX=YO 
DGIRTH=SCRT(0Y*CY*Z1*Z1  ) 

GO  TO  15 
CCNT  INUE 
DY*CY  «Z2/DZ 
DA  =0Y  »Z2/2 • 

YAVG=Y2-DY/2. 

DA YM  = DA«CY2-DY/3.  ) 

CAZM*CA«Z2/3. 

Y0=Y2-DY 

IF (YO.LT.YMIN>YMIN=YO 
DGIATH»S0RT(0V»DY»Z2«Z2) 

15  CONTINUE 
A » A ♦ OA 
B*fi  *0  Y 


V 
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A Y H = A YK*DA YM 
AZN*AZMQAZM 
BYM*3YM*DY»YAVG 

BYI*BYI*0Yn(YAV6»YAVt*0Y»0Y/12«) 
6IRTH-6IRTH+DGIRTH 
14  CO  NT  I NUF 
Y 1 * Y 2 
Z1  =Z2 

10  CONTINUE 

BUP* YMAX -YN  IN 

IF (BWP.LE. 0.0001 )BWP*0. 

AZM=A«ZHIGH-AZM 

RETURN 

ENO 
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Flowchart  of  Subroutine  SECT 
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IF (71 tCT.J.rOOl .ANO.ir ,IE ,0.0001 >60  TO 
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THAXr TO 
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Lilting  of  Subroutine  ENDPT 

SUBROUTINE  f NQP  T ( I j NS  T A«X#GPAFT*XQ«JO) 
OIPFNSICN  DPAFT(25)#X(2!>> 

X * 1 

IF<I.EQ.NSTA)K*-1 
XG*X  ( I ) 

JO*  I 

IF (CPAFTC  I ).GT.O.  )GO  TO  3 
DO  1 J*  2 .»  NS  T A 
JO*I ♦*•< J-1  ) 

IF(CRAFT(JC).GE.O.)GO  TC  2 

1 CONTINUE 

2 CONTINUE 
J1 *JU-K 

X0=(X(J0)>X(J1 ))/2. 

DY*DRAFT (JC)-OPAFT < J 1 ) 

IF (AaSCDY ).LT. 0.0001 >GO  TO  3 
XO  = X (JO)-QMFT (JQ)*(X(JG  )-X(J1  ))/DY 
2 CONTINUE 
PETUfiN 
END 
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Flowchart  of  Subroutine  ENDPT 
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V 


If  <C»  *F  T < JO  ) .fct  .0.  )G0  TO 


Lilting  of  Subroutine  INTEG 

SUBROUTINE  I NTE G ( X t Y » X S T AR T* XfcNDi I START t IENDiR » R1 iR2) 
OINENSION  X(25)»Y(25  ) 

R *0  • 

R1  *0  • 

R2  *0  • 

Y 1 *0  • 

X1-XSTART 

DO  1 I-ISTART# IEND 
DX  *X 1 -X ( I) 

D Y * Y ( I)-Y1 
R*ReDX«(Yl»DY/2.  ) 

R1 -R 1tOX»Yl»(X1-OX/2. )*rx»OY /2.«<X1 -2.»DX/3«) 

R2*R2*0X»«2"Yl/12»^0X»«2»DY/36e*DX»Vl*(X1-DX/2*)»«2*DX«DY/2** 
I (XI -2.«0X/3.)«"2 
X 1 * X < I ) 

Y1 =Y  ( I ) 

1 CONTINUE 

OX*X ( IE  ND  J-XEND 
R*R*0X«Y(IEND>/2. 

ft1*R1eDX«Y(  IfcNO)t*(X(  IEND)-DX/3«)/2. 

R2=R2*DX«»3«Y(  IEN0)/36#*DXiiY  ( IEND)/2.»(X(  IEND  )-DX/3.  )»«2 

RETURN 

END 
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Flowchart  of  Subroutine  INTEG 


• SUBROUTINE  INTEC(X#Y#XSTARTrXEND»  / 

• ISTAAT#  IENDrRrM  2 > / 



I 

I 


DIMENSION  XU5)»Y<25) 

F=a. 

*1=0. 

*2*0. 

Y 1 *0  • 
X1*XSTART 


I 

I 


>:  00  1 1=  ISTAKTr  IEND  s 

<•••••••#•••••••••*••••••••••••••••••••#••••* 

I 

I 


s PX*X1-X<I>  : 

: 0Y*Y<  I ) — Y 1 

: R=R*DX»CYl*DY/2.)  J 

: R1*R1.DX»rl  »<X1-DX/2.)*DX«DY/2.«(Xl-2.»  s 

* DX/3 • ) t 

: K2*R2«-DX*»3»Y1/12.*0X«»*3»DY/36.>DX»'Y1«(  s 

s X1-DX/2.)»«2*DX«*DY/2.*<X1-2.»DX/3.)**2  s 

s X 1 =X ( I ) : 

s Y 1 = Y ( I ) : 

>•••••••■•••••••••••••••••••••••••••••••••••: 

1 


I 1_  I 


: CONTINUE  s 

>•••••••••••••••••••#•••••••••• •••••••••••••* 

I 

I 

: 0X*X( IENO)-XENO  J 

: R*R^GX« Y( IEND )/2.  : 

: R1*M*CX»Y(IEND)»(X(ItN0)-DX/3.)/2.  s 

: R2=R2*PX*»3»Y(IEND )/36.+DX«Y( IFN0)/2.«(X(  : 

: IENO >-GX/3. >«»2  : 

<•••••••••••••••••••••••••••••••••••••••••••: 

I 

I 

/•• i • i •••••••  i ••  i ••  i . 

: EETUAN  : 

■••••••••••••■•••#••/ 
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APPENDIX  B - DEFINITION  OF  PROGRAM  VARIABLES 


The  following  variables  are  used  in  the  main  program.  (Following  these  are 
definitions  of  variables  used  in  the  subroutines.) 


A 

AMAX 

AWP 

AYC 

AYM 

AYM1 

AZC 

AZM 

AZMI 

A1 

B 

BEAM 

BLANK 

BML 

BMT 

BWP 

BWP1 

BYC 

BY  I 

BYI1 

BYM 

BYM1 

B1 

CA1 

CA2 

CB 

CIL 


CIT 

COSPHI 

CP 

CVP 

CW 


CX 


an  array  of  section  areas 
the  maximum  section  area 
the  area  of  the  waterplane 

an  array  containing  the  transverse  center  of  the  section  areas 
an  array  containing  the  transverse  moments  of  the  sections 
a temporary  variable  for  AYM 

an  array  containing  the  vertical  centers  of  the  station  areas 
an  array  containing  the  vertical  moments  of  the  station  areas 
a temporary  variable  for  AZM 
a temporary  variable  for  A 

an  array  containing  the  waterplane  width  of  the  various  sections 
the  maximum  submerged  beam  of  the  boat 
a dummy  variable 

the  longitudinal  metacentric  height  above  the  center  of  buoyancy 
the  transverse  metacentric  height  above  the  center  of  buoyancy 
an  array  containing  the  waterplane  beam  at  each  station 
a temporary  variable  for  BWP 

an  array  containing  the  transverse  center  of  the  waterplane 
at  each  station 

an  array  containing  the  transverse  moment  of  inertia  at  each 
station 

a temporary  variable  for  BYI 

an  array  containing  the  transverse  moment  of  the  waterplane 
at  each  station 
a temporary  variable  for  BYM 
a temporary  variable  for  B 

an  array  containing  a coefficient  defined  as  the  section  area 
divided  by  the  maximum  section  area 

an  array  containing  a coefficient  defined  as  the  section  area 
divided  by  the  area  of  a rectangle  with  the  same  beam  and  draft 
the  block  coefficient  based  on  waterline  length,  beam  and 
maximum  draft 

the  longitudinal  inertia  coefficient  defined  as  the  longitudinal 
waterplane  inertia  divided  by  the  inertia  of  a rectangle  with 
the  same  length  and  beam  as  the  waterplane 
the  transverse  inertia  coefficient  defined  similar  to  CIL 
a temporary  variable  containing  the  cosine  of  the  angle  phi 
the  prismatic  coefficient  defined  as  the  block  coefficient 
divided  by  the  transverse  area  coefficient 
the  vertical  prismatic  coefficient  defined  as  the  block 
coefficient  divided  by  the  waterplane  coefficient 
the  waterplane  coefficient  defined  as  the  waterplane  area 
divided  by  the  area  of  a rectangle  with  the  same  length  and 
beam  as  the  waterplane 

the  transverse  area  coefficient  defined  as  the  maximum  section 
area  divided  by  the  area  of  a rectangle  with  the  same  beam  and 
draft  as  the  maxlmums 
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DISP 

DPHI 

DRAFT 

DRMAX 

DT 

DTHETA 

GIRTH 

GIRTH1 

GML 

GMT 

GZL 

GZT 

I 

I END 
INPUT 

I PRINT 
IPTS 
ISTA 
I START 
ITEST 

II 
12 
J 
K 
K1 
L 

LCB 

LCBM 

LCF 

LCFM 

LCG 

LCG0 

LMOM 

LWL 

LWPI 

NCASES 

NOFFS 

NP 

NPHI 

NPTS 

NS 

NSTA 

NS1 

NT 

NTHETA 

NUMBER 

N1 

N2 

N3 


■ the  displacement  of  the  hull  in  pounds 

• an  incremental  value  of  the  angle  phi 

- an  array  containing  the  draft  at  each  section 

• the  maximum  draft 

■ an  incremental  value  of  the  draft 

■ an  Incremental  value  of  the  angle  theta 

• an  array  containing  the  girth  of  each  station 

• a temporary  variable  for  GIRTH 

- the  longitudinal  metacentric  height  above  the  center  of 
gravity 

• the  transverse  metacentric  height  above  the  center  of  gravity 

- the  longitudinal  righting  arm 

- the  transverse  righting  arm 

- a counter 

- the  subscript  of  the  aftermost  submerged  station 

• a variable  defining  the  device  number  for  input  of  data  cards 

- a variable  defining  the  device  number  for  the  printout 

- a temporary  variable  for  NPTS 

- a counter 

- the  subscript  of  the  forvardmost  submerged  station 

- a variable  used  in  testing  for  equilibrium 

■ a counter 

- a counter 
“ a counter 

- a counter 

- a counter 

- a counter 

- the  longitudinal  center  of  buoyancy 

- the  longitudinal  moment  of  buoyancy 

- the  longitudinal  center  of  flotation 

- the  longitudinal  moment  of  the  waterplane  area 

■ the  longitudinal  center  of  gravity 

- the  input  value  of  the  longitudinal  center  of  gravity 

• the  longitudinal  righting  moment 

■ the  waterplane  length 

- the  longitudinal  waterplane  Inertia 

- the  number  of  cases  to  be  calculated  for 

• the  number  of  offsets  to  be  read  in 

- a temporary  variable  for  NPTS 

“ a control  code  defining  whether  the  heel  angle  should  be 
allowed  to  seek  equilibrium 

■ an  array  containing  the  number  of  points  on  each  station 

- a variable  used  in  setting  up  a table  ^or  output  of  the  results 
« the  number  of  stations 

- a variable  used  in  setting  up  a table  for  the  output  of  results 

■ a control  code  defining  whether  the  draft  should  be  allowed  to 
seek  equilibrium 

- a control  code  defining  whether  the  trim  angle  should  be 
allowed  to  seek  equilibrium 

■ an  array  containing  the  alphanumeric  station  numbers 

- a counter 

• a counter 

■ a counter 
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PHI 
PHID 
PH  10 
RHO 

SINPHI  - 
T 

TCB 

TCBM 

TCB1 

TCF 

TCFM 

TCF1 

TCG 

TCG0 

TCG1 

THETA  - 

THETAD  - 

THETA0  - 

THETA1  - 

TITLE  - 

TMOM 

TWPI 

T0 

T1 

VCB 

VCBM 

VCB1 

VCF1 

VCG 

VCG0 

VCG1 

VDISP  - 
VFORCE  - 
VOLUME  - 
WEIGHT  - 
WS 
X 

XEND 
XSTART  - 
Y 
Y0 
Z 

ZHIGH  - 

ZLOW  - 

Z0 


the  heel  angle  in  radians 
the  heel  angle  In  degrees 
the  input  value  of  heel  angle 

the  density  of  the  water  in  pounds  per  cubic  foot 
the  sine  of  the  heel  angle 

the  depth  of  water  relative  to  the  coordinate  system 

the  transverse  center  of  buoyancy 

the  transverse  moment  of  buoyancy 

a temporary  value  for  TCB 

the  transverse  center  of  flotation 

the  transverse  moment  of  the  waterplane  area 

a temporary  value  for  TCF 

the  transverse  center  of  gravity 

the  input  value  of  the  transverse  center  of  gravity 

a temporary  value  of  TCG 

the  trim  angle  in  radians 

the  trim  angle  in  degrees 

the  input  value  of  trim  angle 

a temporary  value  of  THETA 

an  alphanumeric  title  which  is  used  in  the  output 
the  transverse  righting  moment 
the  transverse  waterplane  inertia 

the  input  value  of  draft  relative  to  the  coordinate  system 

a temporary  value  of  T 

the  vertical  center  of  buoyancy 

the  vertical  moment  of  buoyancy 

a temporary  value  of  VCB 

the  vertical  center  of  flotation  in  original  coordinate  system 
the  vertical  center  of  gravity 

the  input  value  for  the  vertical  center  of  gravity 

a temporary  value  for  the  vertical  center  of  gravity 

the  displaced  volume 

the  residual  vertical  force 

the  desired  displaced  volume  for  equilibrium 

the  weight  of  the  boat  in  pounds 

the  wetted  surface 

an  array  containing  the  longitudinal  offsets  of  the  stations 
the  longitudinal  coordinate  of  the  aftermost  point  of  the 
waterplane 

the  longitudinal  coordinate  of  the  forwardmost  point  on  the 
waterplane 

a two-dimensional  array  containing  the  transverse  coordinates  of 
the  offsets  defining  the  hull  form 

a two-dimensional  array  containing  the  input  values  of  the 
transverse  coordinates  of  the  hull  form 

a two-dimensional  array  containing  the  vertical  coordinates 
of  the  offsets  defining  the  hull  form 

a temporary  value  defining  the  vertical  location  of  the  water 
surface  at  each  station  relative  to  the  coordinate  system 
an  array  containing  the  lowest  points  on  each  station 
a two-dimensional  array  containing  the  input  values  of  the 
vertical  coordinates  defining  the  hull  form 
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The  following  variables  are  used  In  subroutine  SECT.  (Variables  referenced 
as  subroutine  arguments  are  not  Included  In  these  definitions;  they  are  the 
same  as  In  the  main  program.) 


DA 

DAYM 

DAZM 

DGIRTH 

DY 

DZ 

J 

YAVG 

YMAX 

YMIN 

Y0 

Y1 

Y2 

Z1 

Z2 


an  Incremental  value  of  ch«  ut« 

an  incremental  value  of  the  ransverse  moment  of  area 
an  Incremental  value  of  the  vertlial  moment  of  area 
an  Incremental  value  of  the  girth 

the  transverse  distance  from  one  point  to  the  next 
the  vertical  distance  from  one  point  to  the  next 
a counter 

the  average  of  two  transverse  offsets 
the  largest  transverse  offset  encountered 

the  smallest  transverse  offset  encountered,  including  negative 
values 

a temporary  variable  for  a transverse  distance 
a temporary  variable  for  a transverse  offset 

a temporary  variable  for  a transverse  offset 

a temporary  variable  for  a vertical  offset 

a temporary  variable  for  a vertical  offset 


The  following  variables  are  used  in  subroutine  ENDPT.  (Variables  referenced 
as  subroutine  arguments  are  not  Included.) 

DY  ■ an  Incremental  distance 

J -a  counter 

J1  «a  counter 

K "a  value  defining  whether  the  forward  or  after  end  point  of 

the  water  plane  is  being  investigated 


The  following  variables  are  used  in  subroutine  INTEG.  (Variables  referenced 
as  subroutine  arguments  are  not  included.) 

DX  * an  incremental  distance  in  the  X direction 

DY  * an  incremental  distance  in  the  Y direction 

I - a counter 

XI  -a  temporary  variable  for  the  X array 

Y1  -a  temporary  variable  for  the  Y array 
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APPENDIX  C - SAMPLE  RUN 
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* 


VALUES  n FAKEKThESIS  ARE  IN  NON -ROTATED  ORIGINAL  COORDINATES 
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PBS  HYCROSTATICS  PROGRAM  COAST  GUARD  R AND  D CENTER 


• 

• *ESFA*CH 

ft 

BOAT  75-15 

12  FT 

JONBOAT 

ft 

• T ( 

• 2 u6  ) 

• 1 94 

M H « « 

FT 

•DISPLACEMENT 

463.992 

LB 

ft 

•THETA  ( 

.COO ) 

• 30C 

DEG 

•VOLUME  DISP. 

7.436 

CU  .F  T 

• P H I ( 

c u «C  DC  ) 

O.GOC 

DEG 

• LCB 

- • 8 6 C 

FT 

• RE  I GHT 

46  4 • LiC  C 

LB 

•TCb  ( -.822) 

-.869 

FT 

• L Lb 

-1 . 5 2C 

FT 

• VC  B ( .284  ) 

••lil  4 

FT 

•TCO  ( 

.CCD) 

-.51c 

FT 

•watefplane  AREA 

24.1 31 

SQ.FT 

• VCu  ( 

1.510) 

1 .41V 

FT 

«L  CF 

-.653 

FT 

• 

•TCF  ( -.576) 

-.634 

F T 

II 

ft 

» V CF  ( .416) 

« 

• 194 

FT 

• UA  TERPL  A NE 

L E NGTh 

1 1 .374 

KINtt 

FT 

•LONG. MOMENT 

-306.421 

FT  .L  8 

ft 

• W A TtRPL  A NE 

E-tA** 

2 .45 4 

FT 

•LONG .52 

• «6  60 

FT 

•MAXIMUM  DRAFT 

• 6 4 7 

FT 

ft 

• MA x .SECT  I ON  AREA 

.794 

SO. 

FT»TRANS.MOMkNT 

-163.743 

FT. Lb 

•FRO  WATENPLANE  end 

5.513 

FT 

•TRANS. GZ 

-.353 

FT 

•AFT  WATckPLANE  f NO 

-5. 860 

FT 

II 

• WETTED  SUKFACE 
ii 

28.336 

SC. 

FT»VEATICAL  FOACE 

it 

-.00* 

LB 

• LONG.  INEF  I IA 

21 9,856 

OM 

FT4 

•ELOC*  COEFFICIENT 

.41 1 78 

ft 

■ B ML 

29.567 

FT 

• MAX. SECTION  COE  F 

•5JCCQ 

»«ML 

28.134 

FT 

• PRISMATIC  COE  F 

.823 

56 

ft 

• MATE RPLANE  COEF  . 

.36472 

• TRANS. INfcRT  IA 

10.51 1 

F T 4 

• VERT. PRISMATIC  COtF 

.47620 

• BMT 

1 .41  4 

FT 

• LONG.  INER T IA  COEF 

.06090 

» T 

M 

- .Q2C 

FT 

• TMKS#  INERTIA  CO^P 

ft 

.06257 

• VALUES  IN 

PARENTHESIS  ARE  IN  NON-ROTATED  ORIGINAL  COORDINATES 

ft 

ft 
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RBS  HYDROSTATICS  PROGRAM 


COAST  CUAFC  R AND  D CENTER 


• RESEARCH  BOAT  75-15  12  FT  JONBOAT 


■ T ( 

• 1 E2  ) 

.171 

«•*»•» 

FT 

•DISPLACEMENT 

464.718 

LB 

» II 

•theta  < 

-1  .462  ) 

-1 .394 

DEG 

•VOLUME  DISF. 

7.447 

CU.F  r 

• PHI  ( 

2U.CCC  ) 

O.rico 

DEG 

• LCB 

-1  .553 

t T 

■ RE  IGHT 

46 4 • JG  G 

LB 

•TCB  ( -.509) 

-.859 

f T 

« LCD 

-1  .520 

FT 

•VCb  l .239) 

-.00  5 

FT 

• TCb  ( 

.COO) 

-.516 

FT 

•WATERPLAKE  ARtA 

23.55C 

SQ.FT 

• VCG  ( 

1 .510) 

1 .419 

FT 

• LCF 

-1.101 

FT 

II 

• TCF  ( -.57s) 

-.677 

FT 

II 

II 

■ V CF  ( #395 ) 

» 

.171 

FT 

• m A TERPLAM 

LEM  Th 

11.261 

ft  » II 

FT 

• LONG .HOME  NT 

- .61  9 

FT  .L  B 

II 

• RATER  PLANE 

E t AM 

2.71  6 

FT 

•L0N0.G2 

-.002 

FT 

•MAXIHUM  DRAFT 

.76  2 

FT 

H 

• MAX  .SECT  I GN 

AREA 

1 .U35 

SO. 

FT»TRANS. MOMENT 

-159.007 

F T .L  3 

• FRO  bATEPPLANE  END 

5.36  3 

FT 

•TRANS. G2 

-.34  2 

FT 

•AFT  KATERPLANE  END 

-5.65E 

FT 

II 

•WGTTEO  SURFACE 

« 

27.587 

SO. 

FT»\ERTICAL  FO*CE 
« 

.71  8 

LP 

• LONG. INERTIA 

200.366 

IMK 

F T 4 

•ELOCK  COEFFICIENT 

.31 951 

II 

• BML 

26.902 

FT 

•max. Section  coef 

• 50  QOG 

* GML 

25.479 

FT 

•PRISMATIC  colf 

.63901 

« 

• watekplane  coef  . 

.77007 

•TRANS .INERTIA 

1 1 • 25C 

FT4 

•VER T. PRISMATIC  COEF 

.41490 

• BMT 

1 .51  1 

FT 

• LONG • INER  T IA  COEF 

.05167 

• GMT 

• 066 

FT 

•TRANS. INERTIA  COLF 

.04966 

• VALUES  IN  PARENTHESIS  ARE  IN  NON-ROTATED  ORIGINAL  COORDINATES 
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KBS  hydrostatics  program 


COAST  GUARD  R AND  0 CENTER 


RESEARCH  60 A I 75-15  1 2 FT  JONEOAT 


AT  ( .24 4) 

• 2 3t 

ft  ft  ft 

FT 

•DISPLACEMENT 

463.107 

L D 

ft 

•THETA  ( -1.195) 

-1  .163 

DEG 

•VOLUME  OISP. 

7.422 

CU.F  T 

» Prt  I ( -13.244) 

0.000 

OEG 

• L C6 

-1.54  7 

FT 

•WE ISHT 

464 .000 

LB 

• TCB  C 

• 6 C2  ) 

• 63  6 

FT 

• LCg 

-1 . 52C 

FT 

• VCb  ( 

.225) 

.064 

FT 

■ TCG  ( .300 

• 62E 

FT 

•WATERPLANE  AREA 

26.925 

S3.F  T 

• VC»  ( 1.510) 

1 « 4C 1 

FT 

• LCF 

-1 .065 

FT 

ft 

• TCF  ( 

• 2 69) 

.553 

FT 

ft 

ft 

ftVCF  C 

ft 

.312) 

.238 

FT 

• UA  TERPLA  N£  LENGTH 

11 .063 

1 ft  ft  ft 

FT 

• L ONli.MOME  NT 

• 2 9 6 

FT.Lo 

•WATERPLANE  6i  AM 

3.0  69 

FT 

■ L ON  G «GZ 

• J01 

FT 

•MAXIMUM  ORAFT 

.656 

FT 

• MAX.StCT  ION  AkE A 

1 .046 

SO. 

FT.TRANS. MOMENT 

.015 

FT. LB 

• F hi 0 WATERPLANE  END 

5.235 

FT 

•TRANS. G7 

• 000 

FT 

• AFT  KATE  F FLA  N£  END 

-5 • o4£ 

FT 

ft 

■WETTED  SUKFACt 

ft 

30.490 

SO. 

FTft VERTICAL  FORCE 

M 

-.653 

Lb 

•LONG. INERTIA 

223.677 

» ft  ft 

F T 4 

•BLOCK  COEFFICIENT 

.33045 

ft 

• 8M 1. 

30.166 

FT 

•MAX. SECTION  COEF 

.51  636 

• GML 

28.649 

FT 

•FRISMATJC  COEF 

.63992 

ft 

•WATERPLANE  COEF. 

.78667 

■TRANS.INERTIA 

1 7.J40 

FT4 

• VER  T .PR  I SMAT  IC  CCitF 

• 4 £ 0u4 

« BMT 

2.2  96 

FT 

•LONG.INERTIA 

COEF 

.05324 

• GMT 

.975 

FT 

•TRANS. INERTIA 

COEF 

.05218 

» 


VALUES  IN  PARENTHESIS  ARE  IN  NON-ROTATEC  ORIGINAL  COORDINATES 
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APPENDIX  D - TEST  RESULTS 


One  of  several  tests  conducted  for  program  verification  is  presented  here. 
This  test  consisted  of  loading  a jonboat  to  the  load  specified  in  the 
fourth  sample  case  of  Appendix  C,  and  measuring  the  resulting  draft, 
trim,  and  heel.  Trim  and  heel  angles  were  determined  based  on  freeboard 
measurements . * 


BOAT  TEST 

Draft  (feet  amidships)  0.242 
Trim  (degrees)  -1.52° 
Heel  (degrees)  12.45° 


COMPUTER  RUNS 

0.244 

-1.195° 

13.244° 
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